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Abstract 

Top polarization is an important probe of new physics that couples to the top 
sector, and which may be discovered at the 14 TeV LHC. Taking the example of 
the MSSM, we argue that top polarization measurements can put a constraint on 
the soft supersymmetry breaking parameter At- In light of the recent discovery of 
a Higgs-like boson of mass ~ 125 GeV, a large A t is a prediction of many super- 
symmetric models. To this end, we develop a detector level analysis methodology 
for extracting polarization information from hadronic tops using boosted jet sub- 
structure. We show that with 100 fb _1 of data, left and right 600 GeV stops can 
be distinguished to 4c, and 800 GeV stops can be distinguished to 3a. 
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1 Introduction 

1.1 Motivation 

Top physics is an important probe of theories of new physics at the TeV scale, as many 
of these theories posit TeV-scale partners to the top quark in order to resolve the Higgs 
hierarchy problem. These theories in general have chiral structure, thus a measurement 
of the top polarization from the decays of top partners can establish useful constraints 
on them. 

In the case of supersymmetry with .R-parity, the composition of the scalar top partner 
"stop" t in terms of the weak eigenstates and can be constrained by observing 
the polarization of tops in the decay t — > tx±- The fermionic top partners in extra- 
dimensions theories with KK parity and little Higgs theories with T-parity have the 
analogous decays — > tB^ and T" — >■ tB^. As discrete parities are desirable because 
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they limit dangerous contributions to electroweak precision variables pQ and admit WIMP 
dark matter candidates (such as xl, B^ and Bh), the t + $T collider signature provides 
a useful handle on a broad class of well-motivated TeV-scale theories. 

There are also numerous theories with and without top partners containing extra 
massive gauge bosons Z' with decays such as Z' — > ti. Top polarization measurements 
can constrain the chiral structure of their couplings to the top quark. In general, it may 
be produced in new heavy resonances. 

To date there have been a number of studies on measuring top polarization at the 
Large Hadron Collider (LHC). Ref. [2] performed a signal-only Monte Carlo-level anal- 
ysis in the context of gluino decay, and Ref. [3] performed a Monte Carlo-level analysis 
for KK gluons. For the t + class, Ref. pt] performed a signal-only parton-level for- 
mal calculation, followed by Refs. [5] and [6] performing a Monte Carlo-level analysis 
with backgrounds, acceptance cuts and smearing effects. For heavy resonances, Ref. [7] 
considered the P^+jets background and smearing effects, and later Ref. [8] performed a 
signal-only Monte Carlo-level analysis. Ref. [10] elucidated the benefit of jet substruc- 
ture for measuring the polarization of boosted hadronic tops for both classes of theories, 
although without backgrounds and at Monte Carlo-level. 

In this paper we study top polarization for the t + class of theories at detector 
level including all contamination sources (e.g., ISR/FSR and MPI), relevant detector 
effects (e.g., magnetic field) and backgrounds. To this end we focus on pair production 
of 600 GeV and 800 GeV t\ at the 14 TeV LHC under the simplified model in which 
they decay entirely to tXi, where Xi — B and m x o = 100 GeV. We choose to focus on 
supersymmetry also because it may be the most well-motivated of this class of theories, 
solving the hierarchy problem up to Planck scale as well as enhancing gauge coupling 
unification at high scale. 

First we will briefly review the kinematics of top polarization and the phenomenology 
of tL-tn mixing. Then we will describe our simulation and analysis methodology, and 
present our results for the expected sensitivity to the stop mixing angle. Finally, we will 
look ahead to possibilities for improving and ramifying our methodology. 

1.2 Top polarization 

Measurement of top polarization is possible because top quarks undergo weak decay 
prior to hadronization, so the top decay products carry information on the polarization 
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of the parent quark undisturbed by the hadronization process. The kinematics of top 
polarization is presented in Refs. [21 HJ 13 [9] . The decay products of the top quark have 
the angular distributions 

1 dT /T"\ i n /i\ 

rT( fl~T K 1 + Vtk f cos9 tf (1) 

Td(cos9 t f) 

where cos 9 t f is the angle between the daughter momentum and the top spin axis in the 
top rest frame; we can take the latter as the direction of top momentum in the lab frame. 
Vt — ±1 is the polarization of the top quark, and kf is the "spin analyzing power" of the 
daughter flavor. For the 6-quark, 

h = - m \~ m \ c -0.4 (2) 

whereas for the lepton daughter of a leptonic top decay one finds ki — 1. Consequently 
we can measure Vt by observing the distribution of cos8 t f. For the case of hadronic tops, 
this can be done directly by reconstructing the top and resolving the 6-quark daughter. 
However, for leptonic tops, one cannot fully reconstruct the top momentum due to the 
neutrino from the W decay. It has been proposed to define cos Q t i i n an "approximate 
rest frame" in semileptonic events with reconstruction of the accompanying hadronic 
top [5j, or require the leptonic top in a semileptonic event to be highly boosted such 
that one can use alternative measurables which are insensitive to top momentum in the 
limit j3 — > 1 [HE]- Yet another proposal is to use judicious cuts to preserve polarization 
information in the lab frame despite the boost of the leptonic top [8]. 

In our analysis we measure the polarization of hadronic tops, which has not only the 
advantage of being more simple than leptonic methods, but also that of greater statistics, 
as 89% of top pairs have one hadronic top whereas only 44% of top pairs are semilep- 
tonic pj]. Moreover, leptonic top analysis entails all the difficulties of identifying isolated 
leptons in a real hadron collider environment. Standard cone-based lepton isolation lose 
efficiency with increasing boost, hampering polarization measurements for heavy parent 
states. However, this may be ameliorated by narrowing the isolation cone size as lepton 
p T increases [H [T2] . 

1.3 Stop mixing 

As written in the review [15] . the stop mass matrix in the weak basis under the minimal 
supersymmetric standard model (MSSM) is 

£ m{ = -(t* L t* R ) m ?(f) (3) 



m ? = ( m Q 3 + m t + A u L V ( a *t sin P ~ m COS (3) \ 

1 \ v(a t s'm /3 - /j,* y t cos /3) m? 3 + m\ + A SjJ J 

A^ L = Q - sin 2 ^ cos(2/3)m| (5) 
2 

A Sfl = - sin 2 9 W cos(2/3)m| . (6) 

If at and /x are real, the mixing can be represented by the rotation 

t\ \ _ ( cos 9 — sin 9 \ ( 
t 2 J \ sm9 cos9 J \ Ir 

where 

2m t {A t — fi/ tan/3) 



(7) 



tan 29 



m Q 3 _m l 3 + Aa* - A,- 



and in which m t = y f f sin/3 and ^ = at/yt have been substituted. 

Let us consider the case that tan (3 S> 1 and /x, mg 3 and of the same order at 

low scale; we also assume toq 3 > m Ui at low scale due to renormalization. Then if A t = 
we see that t\ ~ (0 ~ 7r/2); conversely, if At > /x then ti will have a significant 
component {6 ~ 7r/4). If indeed the mixing angle # can be measured, then under these 
assumptions A t may be strongly constrained. 

These assumptions can be eased by incorporating other measurables. For example, 
knowledge of the t\ mass (e.g., from production rates or kinematic constraints) and the 
Higgs mass, which takes large corrections at one-loop that depend on A t , /x, tan (3, and 
mi 2 [14], would leave only the relationship between niQ 3 and m S3 as a model-dependent 
quantity. The recent discovery of a Higgs-like boson [151 EE] with a mass suggesting a 
large value of A t for TeV-scale supersymmetric scalars [T7| may already be hinting at 
such a correction. 

It has also been proposed to consider the ratios of the branching fractions of stops 
decaying to neutralinos and charginos in order to constrain stop mixing parameters. 
However, it may require input from both a hadron collider and a linear collider to have 
sufficient information [18]. In either case, this may be another set of observables which 
may be useful in tandem with direct polarization measurements. 

To connect the stop mixing angle with top polarization, we consider the interaction 
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term for the decay in our simplified modefjin which Xi — B 

AC = q 



- 9%X\ 



cos9 ( t L + sm6 Q ) t R 



(9) 



where the quantities in the parentheses are the hypercharges. Then the observed "effec- 
tive" mixing angle from the decay t\ — > t\\ is 



tan ^obs = 4 tan 9 



(10) 



Thus for or simplified model the stop mixing angle is amplified in the top polarization 
mixing angle, increasing the sensitivity to tn and reducing the sensitivity to small 
admixtures. Under the assumptions given before, this therefore reduces the sensitivity to 
small values of A t . 

Putting together the Higgs mass correction and top polarization arising from stop 
mixing, we can see the theoretical sensitivity of A t to the measured top polarization in 
Figure [l] Here, the top polarization is defined as (c| — c 2 R )/(c 2 L + c 2 R ), where cl and cr 
are the coupling strengths (1/6) cos# and (2/3) sin# to the left-handed and right-handed 
tops, respectively. Indeed, sensitivity improves as the tops become more left-handedfl 

2 Simulation and analysis 

The goal of our analysis is to show that top polarization information can be obtained by 
reconstructing boosted tops in a realistic hadron collider environment. First we describe 
the Monte Carlo generation of our data, then the detector simulation, and finally the 
reconstruction of top jets using the physics objects from the simulation. 

2.1 Event generation and detector simulation 

Using Herwig++ 2.5.0 [TH] with all physics effects (hadronization, ISR/FSR, MPI) in- 
cluded, we generated left (sin^ = 0), mixed (tan 6* = 0.25) and right (sin 6* = 1) tit\ 
samples with masses 600 GeV and 800 GeV for the 14 TeV LHC under our simplified 



1 For a discussion with completely general Xi, see e.g. Ref. [S]. 
2 For this parameter scan, tan/3 = 10, fx = 1692 GeV and tua 
assumptions in our discussion. 



1791 GeV, consistent with the 
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Figure 1: Constraint on trilinear coupling A t vs. top polarization measurement from the 
MSSM Higgs mass correction for 800 GeV light stop and 123.8 < rrih < 126.2. 



model. We calculated NLO production cross sections using Prospino2.1 [20], shown in 
Table ffl 

We consider the backgrounds tt+jets, Z+jets, VT+jets and tt+Z, which we generated 
using MadEvent/MadGraph 5.1.3 [HI [22] + PYTHIA 6.4.25 [23] (also with all physics 
effects), taking their leading order cross sections which are also shown in Table [I] All 
extra jets are five-flavor (g + u, d, c, s, b). 



Process 


Generator-level cut 


Cross section 


tit*, m = 600 GeV 
trf*, m = 800 GeV 




218 fb 
36.8 fb 


tt+ < 2 jets 
[Z — > up)+ < 3 jets 
(W -> [Z,r]i/)+ < 3 jets 
ti+(Z ->■ vv) 


p T ,ji > 300 GeV 
p T>jl ', $ T > 250 GeV 
p T>jl , $ T > 300 GeV 


40.6 pb 
7.8 pb 
4.57 pb 
0.11 pb 



Table 1: Signal and background cross-sections at 14 TeV LHC. 



For our detector simulation we used Delphes 2.0.2 |24j. We modified the Delphes 
codebase to use Fast Jet 3.0.3 [2S] instead of the bundled version, as the newer version 
has an interface to manipulate sub jets at specific clustering scales or steps. This allows us 
to "prune" the clustering tree to remove contamination, then store the resulting subjets, 
all within the Delphes analysis pipeline. 
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The Delphes detector settings are tuned to ATLAS, with the hadronic calorimeter 
grid set to match that in the ATLAS TDR (26]. The magnetic field is turned on in the 
simulation. 



2.2 Cuts 



We implement the following event cuts, designed to increase the significance S/ V S + B 
for our characteristic t + $t signature: 

1. $ T >300 GeV. 

2. Leading fat jet pt >400 GeV. 

3. If there are no leptons w/ pr > 5 GeV, require subleading fat jet pr >100 GeV. 
This cut suppresses processes like W/Z+jets in which the second fat jet is from 
QCD, as these jets are likely to be soft. Since the majority of signal events have 
leptons due to W — > Iv and the decay of 6-flavored mesons, we require that there 
are no leptons for this cut. Thus, this cut is most effective against Z+jets. 

4. Lepton is not collimated with ~$t- For every lepton with px > 5 GeV, require 



This selects against high ft? arising from boosted leptonic W decays in tt+jets 



smaller in these processes than from a top partner decay. 

5. Hard subjet is not collimated with $t- For every subjet with px > 50 GeV, require 
the same as above. This works against hadronic r from W decays in tt+jets and 
W+jets, as well as highly collimated 6-subjets from top decays in tt+jets. 

6. Require at least one top tagged jet, using the procedure described in the next section. 

7. 225 < M T2 < 650 for 600 GeV stops, and 325 < M T2 < 850 for 800 GeV stops. 
M T2 is calculated for the leading and subleading jet, with m x = 0. For the leading 
jet we used the reconstructed top jet if top tagged, otherwise we used the trimmed 
jet; similarly for the subleading jet. We employed the Mt 2 code of Ref. [37]. 

The resulting cut flow for signal events is shown in Table |2j and for background events 
in Table [3] for 14 TeV LHC @ 100 fb _1 . The cuts reveal some preference for right stops, 
which is noted in Ref. [5]. 




< 0.4 



(11) 



(#r + Pr,{)/(350 GeV) 



and W+jets, as the opening angle between the lepton and $t is likely to be much 



S 



Cut 


Stop 600 GeV 


Stop 800 GeV 


Left 


Mixed 


Right 


Left 


Mixed 


Right 


# 


Pre-cut 


21900 


21900 


21900 


3680 


3680 


3680 


1 


$ T > 300 GeV 


9513 


9739 


9857 


2394 


2411 


2433 


2 


p Td i > 400 GeV 


5415 


5496 


5472 


1816 


1835 


1825 


3 


If m = 0, > 100 GeV 


5150 


5220 


5192 


1756 


1773 


1764 


4 


lepton/^ collimation 


4155 


4315 


4394 


1515 


1555 


1559 


5 


subjet/^-j 1 collimation 


2914 


3046 


3084 


1171 


1209 


1208 


6 


# top tag > 1 


1014 


1065 


1082 


450 


456 


463 


7a 


225 < M T2 < 650 


908 


954 


969 








7b 


325 < M T2 < 850 








364 


369 


374 



Table 2: Cut flow for signal events at 14 TeV LHC @ 100 fb" 1 . 



Cut 


tt+jets 


Z+jets 


W+jets 


tt + Z 


# 


Generator-level 


4.06 x 10 6 


7.8 x 10 5 


4.57 x 10 5 


11000 


1 


$ T > 300 GeV 


1.30 x 10 5 


4.96 x 10 5 


4.36 x 10 5 


815 


2 


prji > 400 GeV 


90503 


2.28 x 10 5 


3.20 x 10 5 


351 


3 


If m = 0, ptj2 > 100 GeV 


88133 


1.01 x 10 5 


2.68 x 10 5 


326 


4 


lepton/^T collimation 


21518 


98441 


69865 


237 


5 


subjet/^r collimation 


4412 


60860 


37852 


149 


6 


# top tag > 1 


1140 


305 


99 


52 


7a 


225 < M T2 < 650 


554 


222 


56 


43 


7b 


325 < M T2 < 850 


275 


141 


45 


30 



Table 3: Cut flow for background events at 14 TeV LHC © 100 fb _1 . 



2.3 Top jet reconstruction 

Many aspects of this analysis are well- reviewed in Ref. [27] . 

2.3.1 Jet clustering and grooming 

Hadrons are clustered as fat jets using the Cambridge/Aachen algorithm [28] with cone 
size R — 1.2 and sub jet cone size AR = 0.2. These numbers are chosen such that the 
distribution in the number of subjets per fat jet peaks at ~ 3 in our signal samples after 
the following grooming procedures: 

1. The jet clustering trees are "pruned" [29] using the mass-drop condition [30] 

m i „_ 1 < 0.8 x m jn (12) 

where rrij n _ 1 is the invariant mass of the hardest parent jet, and rrij n is the invariant 
mass of the child jet at clustering step n. We also require the subjet separation 



9 



condition 

4 T (jn-i,i, jn-1,2) > (Ai?) 2 ■ m 2 jn (13) 

where c4 T is the kx distance between the two parents jets. This removes contami- 
nation from MPI and ISR, improving top reconstruction quality. We take the fat 
jet (and its subjets) at the clustering step where both conditions are satisfied. 

2. This jet is then "trimmed" |31j . removing subjets with p? < 10 GeV. This further 
reduces contamination. 

The trimmed jet is then fed to our top reconstruction algorithm. 

2.3.2 Reconstruction and tagging 

The following algorithm is attempted for every fat jet: 

1. Require pr > 400 GeV for the untrimmed jet. 

2. Require at least three subjets. 

3. Require that one 6-subjet jb is reconstructed in the jet. 

4. Find the subjet combination jxjijbi of which no pair of subjets are within AR of 
each other, and which gives the closest invariant mass to m t . Require also that this 
invariant mass be in the window (150, 200) GeV. 

5. If there is no successful tag, require at least four subjets and retry the step above 
with four-subjet combinations j\j2]z3b- 

6. Optionally require one two-subjet combination to have an invariant mass in the 
loose W mass window (50, 110) GeV. We present our main results both with and 
without this requirement. 

This algorithm differs from the Johns Hopkins top tagger [32] by declustering more than 
two steps in the pruning stage if necessary, using a mass drop condition rather than a px 
drop condition, having an absolute rather than fractional trimming threshold, requiring 
a 6-tag instead of imposing a W mass condition, and not imposing a top helicity angle 
condition (as this is our observable). The algorithm also differs from the CMS tagger [33] 
by not requiring a minimum two-subjet invariant mass. 
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Kinematic region 


Efficiency 


p T < 30 GeV 


0% 


30 GeV < vt < 60 GeV 


60% 


60 GeV <p T < 450 GeV 


70% 


450 GeV <p T < 1000 GeV 


60% 


1000 GeV < p T 


0% 



Table 4: Utilized fe-tagging efficiencies. 



Our algorithm differs also from the HEPTop Tagger [31] by not imposing the various 
two-subjet mass requirements, having a mandatory 6-tag, and also by implementing four- 
subjet reconstruction. 

We make these choices to enhance top tagging efficiency, presuming that the top parent 
particle has already been discovered. Nonetheless, top mistagging does not overwhelm 
the signal as will be apparent in our results. 

2.3.3 6-tagging subjets 

We require 6-tagging in our top reconstruction algorithm to reconstruct the observable 
cos 9 t b with high fidelity, but also to suppress top mistagging from background processes. 
Utilizing recent advances in 6-tagging by the LHC detector collaborations, we employ the 
6-tagging efficiencies (shown in Table [4]) recently validated at 7 TeV LHC by CMS for 
their CSVM tagger |35J. We impose the upper limit px = 1000 GeV to be conservative, 
though it is not indicated by the CMS study. We choose to use the CMS efficiencies since 
they are validated up to px = 670 GeV, though a recent ATLAS study [36] shows similar 
efficiencies up to 200 GeV. 

However, we do not implement mistagging, as this depends on various factors that 
are best implemented by experimenters — mistagging rates vary rapidly with tagging 
efficiency and so are sensitive to systematic uncertainties that we cannot model. We 
expect charm jets to contribute to most of the mistags. For example, in Ref. [35] a 60% 
6-tagging efficiency nominally results in a 10% charm mistagging rate in Monte Carlo, 
whereas for 50% 6-tagging efficiency this drops to 4%. The mistagging rate for lighter 
flavors is ~ 1%, so this contribution is negligible. 

We apply these 6-tagging efficiencies at parton level. To reconstruct a 6-subjet we 
sum all the subjets within AR = 0.2 of the 6-parton, as this captures some hard FSR. 
If more than one 6-parton yields a matching subjet inside a given fat jet, that fat jet is 
rejected and is not top tagged. 
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3 Results 



3.1 Reconstruction quality 

To illustrate the effectiveness of our reconstruction method, we show the jet invariant mass 
distribution at different stages in Figure [2j For signal processes and tt+jets, few jets are 
lost in the pruning phase; however, without jet cluster pruning, we find that the cos 6tb 
distribution loses fidelity compared to the parton-level expectation. The subsequent 
trimming step is essential for removing soft contamination, resulting in a shift of the 
invariant mass lower towards the correct top mass. Then, requiring at least three subjets 
and a 6-tag narrows the distribution further. Finally, top reconstruction assembles the 
subjets with correct invariant mass. Modulo 6-tagging efficiency, the top tagging efficiency 
for hadronic tops in our samples which pass the px cut is ~ 60%. Conversely, top 
reconstruction suppresses Z+jets and H^+jets by a factor of 0(100). 
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Figure 2: Jet invariant mass sequence for 600 GeV mixed stops (upper left), 800 GeV 
mixed stops (upper right), tt+jets (lower left) and Z+jets (lower right). In each panel, 
from back to front, shown is the jet mass distribution after p? > 400 GeV cut (blue), 
requiring successful pruning (yellow), trimming (gray), requiring three or more subjets 
(cyan), 6-tagging (magenta), and finally top tagging (green). 
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The quality of reconstruction is also apparent in the signal-only cos 6^ distributions 
shown in Figure [3] One sees that the parton-level and reconstructed distributions coincide 
up to statistical fluctuations. The distributions at pre-cut/tag parton-level match Ref. jl], 
except they have a downward tilt due to the 6-partons losing energy to FSR. One sees 
also that there is less contrast between left and mixed stops with mass 600 GeV than 
with mass 800 GeV. This is due to the W daughter subjets not passing the trimming 
threshold. 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.1 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.I 



1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 3: cos Q t b distributions for stops only, normalized to one. Left to right are left, 
mixed, and right stop samples; upper row is for 600 GeV stops, the lower for 800 GeV 
stops. Shown in each panel are the reconstructed distribution (black), parton-level dis- 
tribution (blue), and parton-level distribution before cuts and tagging (green). 



3.2 tt + jets control region 

As experimenters are likely to rely on data rather than Monte Carlo for the tt+jets 
background, we demonstrate a control region by inverting cuts #4 and #5, requiring 
that there be at least one lepton or hard subjet highly collimated with Ifi?. An array of 
cos Q t b distributions for this region is shown in Figure |4j They evince little contamination 
from signal and other backgrounds. 
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It may be of concern that the collimation cuts $4 and #5 might introduce distortions 
to the cos 6tb distributions obtained from top tagging. In Figure [5] we show that the polar- 
ization distribution for tt+jets is invariant under these cuts up to statistical fluctuations. 
We show also in Figure [6] that the My 2 cut does not distort the polarization distribution 
of tt+jets in the control region. Thus, it is safe to use a data-driven tt+jets background 
from this region. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
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-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Figure 4: Event distributions of cos Q t b for tt+jets in the control region at 14 TeV LHC @ 
100 fb _1 . Left-to-right are left, mixed and right stop samples; upper row is for 600 GeV 
stops, the lower row for 800 GeV stops. Shown in each panel top to bottom is signal 
(red), tt+jets (blue), Z+jets (green), I^+jets (yellow), and tt + Z (gray; hardly visible). 



3.3 Sensitivity to stop mixing 

Our measurable is cos 9 t b, where Q t b is the angle between the b daughter and the top spin 
axis in the top rest frame. To calculate sensitivity, we sum the cos Qtb distributions from 
signal and background processes. For the it + jets contribution, we use the control region 
distribution normalized to the total of the signal region distribution. This sum is shown 
in Figure [7j where we have rebinned the data to use only two bins in order to minimize the 
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Figure 5: Distributions of cos 9 t b for tt+jets, normalized to one. Shown are the distribu- 
tion without collimation cuts (black), with the lepton collimation cut only (blue), with 
the subjet collimation cut only (red), and with both collimation cuts applied (cyan). 

trials penalty. The corresponding p-values for distinguishing stop mixing hypotheses are 
shown in Table |5j For 600 GeV stops at 14 TeV LHC @ 100 fb _1 , left and right mixtures 
can be distinguished to better than 4a, and left /right can be distinguished from the 
mixed state to better than 1.5a. For 800 GeV stops, left and right can be distinguished 
to nearly 3a. 



Truth 


Hypothesis (m^ = 600 GeV) 


Hypothesis (m^ =800 GeV) 


left 


mixed 


right 


left 


mixed 


right 


left 

mixed 

right 


1 

0.058 
8.1 x 10" 6 


0.047 
1 

0.032 


2.6 x 10~ 6 
0.030 
1 


1 

0.17 
0.0018 


0.16 
1 

0.25 


0.0014 
0.24 
1 



Table 5: p- values for distinguishing stop mixing hypotheses at 14 TeV LHC @ 100 fb 1 , 
with no W mass condition. 



Except for ti + Z, most of the tops from the background processes surviving the cuts 
and reconstruction are due to extra QCD jets carrying 6-fiavor. The results therefore 
are improved somewhat by imposing a loose W mass condition in the top reconstruction 
algorithm. If we require the mass of some two-subjet combination (not including the 
6-subjet) to be in the window (50, 110) GeV, we obtain the sum plot shown in Figure [8] 
and the corresponding p-values in Table |6j For 600 GeV stops, left and right can now 
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Figure 6: Distributions of cos 9 t b for the tt+jets control sample, normalized to one. Shown 
are the distribution for M T2 < 225 GeV (black), and for M T2 > 225 GeV. 



be distinguished to better than 4.5a, and left/right can be distinguished from mixed to 
roughly 2a. For 800 GeV stops, left and right can now be distinguished to better than 
3a. However, imposing this W mass condition may increase systematic errors from the 
W/Z +jets backgrounds due to the greater variation between the negative and positive 
bins. This occurs because the subjets of the fake Ws have lower pr at larger cos 9 t b, and 
are less likely to satisfy the W mass condition. 

We conclude that stop mixing hypotheses can be distinguished at the 14 TeV LHC 
with ~ 100 fb" 1 of data. 



Truth 


Hypothesis (m^ = 600 GeV) 


Hypothesis (m^ =800 GeV) 


left 


mixed 


right 


left 


mixed 


right 


left 

mixed 

right 


1 

0.025 
9.0 x 10~ 7 


0.018 
1 

0.018 


1.7 x 10~ 6 
0.017 
1 


1 

0.17 

7 x 10- 4 


0.17 
1 

0.15 


5.6 x 10" 4 
0.14 
1 



Table 6: p- values for distinguishing stop mixing hypotheses at 14 TeV LHC @ 100 fb l , 
with a W mass condition. 
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Figure 7: Sum of cos 9tb distributions from different processes for 14 TeV LHC @ 100 fb _1 , 
with no W mass condition. Left to right are left, mixed and right stop samples; upper 
row is for 600 GeV stops, and the lower is for 800 GeV stops. Shown are the contributions 
from signal (red), tt+jets (blue), Z+jets (green), W+jets (yellow), ti + Z (gray). 
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Figure 8: Sum of cos 9tb distributions from different processes for 14 TeV LHC @ 100 fb _1 
with a W mass condition. Left to right are left, mixed and right stop samples; upper row 
is for 600 GeV stops, and the lower is for 800 GeV stops. Shown are contributions from 
signal (red), tt+jets (blue), Z+jets (green), W+jets (yellow), ti + Z (gray). 
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4 Conclusion and outlook 



In summary, after reviewing the motivation for top polarization measurement, the kine- 
matics of top polarization, and the phenomenology of stop mixing, we described a sim- 
ulation and analysis methodology for the t + $t collider signature that can distinguish 
stop mixing hypotheses at 14 TeV LHC with ~ 100 fb _1 of data. 

There are several possible improvements to the methodology. The first and foremost 
is to include polarization information from leptonic decays, perhaps using the techniques 
of [3 H2] • Also if, in consultation with experimenters, the trimming threshold of px = 
10 GeV or the subjet cone size AR = 0.2 can be reduced, this would enhance the 
performance of the top reconstruction algorithm. Finally, it may be useful to implement 
the top reconstruction technique of Ref. [TU] to supplement 6-tagging, especially at large 
boosts. 

Other possible improvements are to: 

• Implement charm mistagging in the 6-tagging algorithm, though this may be non- 
trivial as discussed in Section 12.3.31 

• Use spin-correlated backgrounds, e.g. from ALPGEN [38]. We did produce a small 
sample of ti+jets in ALPGEN but found no difference in the cos Q t b distribution from 
the MadGraph/MadEvent+PYTHIA events used in this analysis. 

• Calculate backgrounds at NLO using, e.g., MC@NL0 [39J to increase accuracy. How- 
ever, this is likely beyond our computational capabilities unless it is used only to 
normalize the leading-order backgrounds. We estimate that the .fT-factor for the ti 
background to be less than 1.2 given our ^ cut [40J. 

In conclusion, the LHC running at 14 TeV may provide new discoveries in the top 
sector, in which case top polarization will be an important tool for constraining this new 
physics. The methodology in this paper may be useful for this purpose. 
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